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Abstract. Optimal transport from the volume measure to a convex combi- 
nation of Dirac measures yields a tessellation of a Riemannian manifold into 
pieces of arbitrary relative size. This tessellation is studied for the cost func- 
tions Cp{z,y) = ^d^\z,y) and 1 < p < cxd. Geometric descriptions of the 
tessellations for all p is obtained for compact subsets of the Euclidean space 
and the sphere. For p = 1 this approach yields Laguerre tessellations for all 
compact Riemannian manifolds. 



1. Introduction 

Given two probability measures fj, and v on some Riemannian manifold M, a 
transportation map from fi to v is a map T : M ^ M pushing /i forward v, i.e. 
T^/i = v. In the theory of optimal transportation one is interested in a special 
transportation map, namely the one minimizing 

c{x, T(x))dii{x) 

M 

the transportation cost among all possible transportation maps, for a given cost 
function c : M x M ^ R. This minimizer, if it exists, is called optimal transporta- 
tion map. 

A voronoi tessellation of a Riemannian manifold M with respect to a discrete set 
of points 5* is a decomposition of M into different cells with center s G S. Each cell 
consists of those points which are closer to its center than to any other center. A 
rather natural extension are the Laguerre tessellations where different centers can 
have different weights, see jLZOSj . 

The theory of optimal transport is a very natural setting for studying such tessel- 
lations. Just take fi as the normalized volume measure, v as a convex combination 
of Dirac-measures and as the cost function c{x,y) = d{x,y), the geodesic distance. 
In this framework the generalization of Laguerre tessellations to other cost func- 
tions like c{x,y) = -dP{x,y) is straightforward (at least to write the tessellation 
down) . 



Some of the phenomenons studied in this paper already appeared in the work of 
Ma, Trudinger and Wang, Loeper, Kim and McCann (see |MT W05[ iLoeOOl IKM07) ) . 
where they studied the regularity of the optimal transport map. In the case above 
for p = 2 and the manifold M being the hyperbolic disc, the resulting "cells" do 
not have to be connected any more. If one now smears the Dirac points slightly 
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one ends up with a measure which is absolutely continuous to the volume measure 
(the density can even be very nice) but with a discontinuous transportation map! 
Therefore, even the transport between two measures being absolutely continuous 
to the volume measure can be rather irregular. 

1.1. General setting. In I McCOl] McCann showed that on a compact connected 
complete Riemannian manifold {M, g) without boundary there is always an optimal 
transport map T pushing some measure /x, absolutely continuous to the volume 
measure vol, forward to an arbitrary Borel measure v on M under the condition 
that the cost function c is superdifferentiable. The map T is the exponential of the 
gradient of a c-convex function $. A function is c-convex iff ($'^)'^ = <i>, where 

= - inf [c(a;,t/) 

yeM 

We consider the special case of /i(-) = l/vol{M)vol{-), the normalized volume mea- 
sure, and v{-) = X]i>o •^i'^a^i ('); with \i '>Q,^\i — 1 and Xi G M for all i, a 
(possibly infinite) convex combination of Dirac measures, and the cost function 
Cp(x, y) — ^dP{x, y) for 1 < p < cx), where d(x, y) denotes the geodesic distance be- 
tween X and y. Then, the optimal transport map can be written down explicitly, i.e. 
the c-convex function $ can be written down explicitly (the apparent p-dependence 
is suppressed in the notation for the optimal transport map): 

dP 

(1.1) ^{z) = max $i(z), with ^^(z) = (xi,z) + 6i, 6i £ M Vz. 

l<i<n p 

In particular, this shows that there is always a solution to the (even very much 
generalized) tessellation problem. 

We are interested in the geometry of the mass being allocated to a certain Dirac 
point. More specific, how does the set T^^{xi), the cell allocated to Xi, look like? 
Does is have to be connected? 

A nice way to look at these sets is to interprete them as intersections of "halfspaces" , 
i.e.: 

T-\x,) = {geAf:<I>,(g)>$,(g)VjVi} 

= f]{qeM:Mq)>'i'j{Q)} 

Thus, a good way to study the geometry of the cells T^^{xi) is to look for proper- 
ties invariant under sections of these "halfspaces" . In this way we can reduce the 
problem from a many point problem to a two point problem. 

This also directly settles the question of the regularity of the boundary of the 
cells T^^{xi) in the case of finitely many Dirac points. To study this regularity, 
it suffices to study the regularity of the boundary of the halfspaces dH!j = {p G 
M : ^i{q) — ^j{q)}. This question is answered by an implicit function theorem 
for subdifferentiable functions, e.g. like in the second appendix to chapter 10 in 
|Vil09j . and yields 

Lemma 1. is (locally) a (n — 1)- dimensional Lipschitz graph. For all q G dHj 
with q ^ {cutlocus{xi)} U {cutlocus{xj)} there is a neighbourhood U of q such that 
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Figure 1. p = 3/2, a = 5 

U n dHj is the (n — 1) dimensional graph of a function which is as smooth as the 
cost function. 

Thus, the boundary of the set T~^{xi) is, as an intersection of sets with locahy 
Lipschitz boundary, itself a set with locaUy Lipschitz boundary up to the points 
of intersection with the different "halfspaces" . In particular, dHj is a /i-nuU set. 
Considering just finitely many Dirac points, we can thus define the cells T~^{xi) to 
be open. In the case of inifinitely many Dirac points, things become less regular. 
For example, take a dense countable subset of the unit disc as Dirac points. As the 
different cells are starlike, see section^ they become stars consisting of many lines. 

We now turn to the question of the geometry of the cells, and especially the 
connectedness. First, we will study the case of compact subsets of M". Here we 
distinguish between convex and non convex subsets because they behave drastically 
different. Then, in section[3]we consider the problem on a curved space, the sphere. 
It turns out that this case is the nicest. We can describe the geometry of the cells 
well and moreover, opposed to the Euclidean setting, the cells are connected. In 
section 2] we deal with the case p — 1 which can be treated for all Riemannian 
manifolds yielding Laguerre tessellations. 

2. Euclidean case 

2.1. Convex subsets of R". In this section, we will look at compact simply con- 
nected subsets fl of R". Let us assume untill explicitly said otherwise that H. is 
convex. In the case p — 2 and finitely many Dirac points it was shown in section 
4 of [Stu09| and also in [AHA92j that the cells T~^{xi) are convex polytopes and 
therefore simply connected. This can be easily seen from the characterization of 
this set as an intersection of halfspaces. Indeed, dHj is just a straight line. Thus, 
both and Hf are convex polytopes. As intersections of convex sets are convex 
sets, Clj^i Hj is still a convex polytope and we are done. However, this is the 
only case where we actually have connected sets. For all other cases of p, the sets 
T^^{xi) do not have to be connected! Indeed, we have the following 




Figure 2. p = 5, a = 100 



Theorem 1. Letn > 2, f7 C K" he convex compact and simply connected. Consider 
the optimal transport problem for fi and v = \5x^ + (1 — A)^^;^ and cost function 
Cp{x,y) with p £ (1,(X))\{2}. Then, H2 is in general not connected. 

Proof. As the measure is a convex combination of just two Dirac- measures, we 
can safely work in M^. We will take xi ~ (1,0) and X2 = (—1,0). In order 
to characterize the set H2 one needs to study its boundary 9i?2- However, the 
boundary is just some level set a of the function 

m{x, y,p) = \xi - {x, y)\P - \x2 - {x, y)\P = {{x + if - y^'^ - {{x - if + y^'^ 

where {x,y) denotes some point in R^. Figure 1 and figure 2 show how these level 
sets look like for certain values oip and a. The figures show the typical behaviour 
of m for p < 2 and for p > 2. They also show, why theorem [1] holds. 

In the case of p < 2 (figure 1) we can draw a straight line from a point in Hf, 
that is on the left hand side of the level set, to a point in H2 on the right hand side 
of the level set which intersects the level sets three times. Completing this line to a 
convex polygon Q. such that Xi and X2 are inside the polygon, we find a convex set 

such that Hi is not connected. To show that this is indeed possible it suffices to 
show that the graph of the level set has a minimum in y = and a change in sign 
in its second derivative for some y > 0. 

The case p > 2 is slightly easier. The graph of the level set has a maximum at 
y = 0. For a large enough, xi G H^. Thus, we can take a segment of the hyper- 
plane X = c for some c > 1 such that (c, 0) G Hf, which intersects with the graph 
of the level set two times, and complete this segment to a polygone analogously to 
the case above. For this, it suffices to show that the graph of the level set has a 
maximum in y = 0. 

By the implicit function theorem, if ^^{xq, yo) ^ then there is a function r]{y) 
such that in a neighbourhood of (a;o,?/o) we have m{r]{y),y,p) = a. This condition 
holds for all x > 1, the case we are interested in, and 1 < p < 00. Moreover, we 
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have that 



df] dm/dy 
dy dm/dx 



{v{y),y) 



and for the second derivative of rj we get 



(2-1) ^ 



dhi _ 1 
dy'^ dm/dx 



dx^ 



dr] 
dy 



^ +2 



d'^m drj d'^rr 
dydx dy dy"^ 



This means, that even though we do not know how ry cxpUcitly looks hke, we do 
know what its derivative at any point in with x > 1 is. And this is enough for 
our purpose. We can check, if rj has an extremum at t/ = and if we know that 
for any fixed a; > 1 the second derivative of r] changes its sign, we are done. Let us 
first check the extremal points of r]. 

Calculating the first derivative of yields: 



py{{-l + xf + y'^) 



-1+1 



■pj/((l + a;)2+y2) 



2^.,2\-l+f 



-p{-l + x) ((-1 + xf + +p{i + x) ((1 + xf + 2/2)-^+* 



2 . „2^-l+l 



which is zero at y = 0. The second derivative at the point [x, 0) is: 



{{-l + xf) 



((1 



-1+1 



-(-1 + x) ((-1 + + (1 + x) ((1 + :r)2)-^+5 



As wc consider the case x > \,p > 1 the denominator is strictly positive. As the 
numerator is negative for p > 2 and positive for p < 2, rj has at the point {x, 0) a 
maximum for p > 2 and a minimum for p < 2. 

This settles the case for p > 2. For p < 2 it remains to show that the second 
derivative of rj has a change in sign for fixed x. The second derivative at the point 
{x, y) is 



6 



MARTIN HUESMANN 



-p{-l + X) ((-1 + X)^ + y2)i(-2+P) + ((1 + + y2)i(-2+p) 



i-2 + p)ym-l + xf+y^) 



,2\5(-4+p) 



+ {i-l + xf+y') 



2\5(-2+P) 



- (-2+p)2/2((l + a;)2+y2) 



2\3(-2+p) 



+ 



? 2 

p y 



(p{-l + X) ((-1 + X)2 + y2)^(-2+p) _ ^ ^) ^ 3^)2 + ^2)^(^2+^)^ 

X ( ((-1 + xf + y^r^-'^"^ - ((1 + xf + y^r^-'^^^y X 

X {{-2+p){-l + xf ((-1 + xf + y^) + ((-1 + xf + y') -^-'^'^ 

- (-2 + + xf ((1 + xf + y^) - ((1 + xf + y') ^^-'^^^ } 

^ 2(-2+p)pt/^ ^ 

+ a;) ((-1 + xf + y2)l(-2+P) + p(l + a.) ((1 + xf + t/2)i(-2+P) 
X (-(-1 + X) ((-1 + xf + y2) + (1 + ^) ((1 + ^)2 + ^2^ i(-4+p)^ ^ 

X (- ((-1 + xf + yr'-'^^' + ((1 + xf + yr'-'^^') ] 

For .X > 1 the first fraction is positive. Thus, we can concentrate on the sum 
of the three terms inside the square brackets. The common denominator of these 
terms is a square. As we are interested in the sign, it suffices to look at the numer- 
ator after the expansion to the common denominator. This yields three terms, say 
u, V and w. As we have already seen that the second derivative of rj is positive at 
{x, 0) (for 1 < p < 2), it suffices to show that the numerator converges to zero from 
below as y tends to infinity to be able to conclude with the mean value theorem 
that there is a change in sign (a null with change of sign). 

The first term reads 

u = [(-2+p)y2((-l + x)2+y2)^(-^+^)-(-2+py((l+x)2+y2)i(-4+P) 

+ ((-l + x)2+y2)^(-2+P)_((i + ,,,)2^,^2)^(-2+p)- ^ 

X (pi-l + x){{-l + xf+ y') ^^-'^'^ - p{l + x){{l+ xf + y') ^("2+^)) ' 

By using functions like g{z) = {x + 1)^, we find the asymptotic u ~ 8p^(p — 
2) (3 — p)xy^P~^. The second term reads 



2 2 

V = p y 



{ii-i+xf+yr'-'^''-{i^+-r+yr'''^''' 

ip -2){x- If {{X - If + y^y^'-'^ -{P-2){1+ xf ((1 + xf + y')^^'-'^ 

+ ((-1 + xf + y') - ((1 + xf + y2) ^(-2+^') ■ 
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Figure 3. Example [T] 
which behaves Hke v ^ p^{p — 2)^x^24:y^P^^'^ . The last term is 

X + X) ((-1 + + ^ ^ ^) ^ ^^)2 ^ 2^2) 

This behaves asymptotically like w ^ I6p^{p — 2)'^xy^'P^^ . In total, we have an 
asymptotic for y tending to infinity like 

u + v + wr^8{p- l){p - 2)p^xy^P-^ - 

as we assumed 1 < p < 2. Thus, there is a change in sign of the second derivative 
of T] for any fixed x > 1 proving the theorem. □ 

2.2. Non-convex subsets of M". If we now remove the restriction of ft being 
convex we are even worse off. Then, it can happen that T~^{xi) has even infinitely 
many components. 

Example 1. Take two points in M.'^, say Xi and X2 from above, and for simplicity 
c{x,y) — ^{x,y) (this also works for other p > Ij. Arrange the constants bi such 
that T{xi) = X2- By induction, we now build a set fl such that T^^{xi) has infintely 
many components. The basic idea is as follows. dH\ is a straight line. Everything 
on the left is transported to X2 everything on the right to xi . Just cut out a set with 
infintely many components on the right hand side of 9-^2 • 

Cut a nice set f2o out of such that some segment of the straight line dH^ is a 
connected subset of the boundary of flo, say G — dH^ H dflo- Then, take a nice 
subset Wi of Hi such that dWi D dH^ =giCG. Set f^i = 17o U M^i . 
Let W2 be a version ofWi scaled down by some factor k/2 and translated such that 
g2 — dW2 n dH^ C G and i'DfxeW2-veWi d(x, y) > k/2. Then set 0,2 = fliU W2. 
Let be constructed and let Wn-^-i be a version of Wn scaled down by the factor 
k/2 and translated such that gn+i — dWn+i ndH^ C G, TJ^ixew„+i.yew„ d(x,y) > 
(k/2)" and mi^^w^+uyaw, d{x,y) > Ylk=j 2(k/2)'=. Then set fln+i = ri„ U Wn+i. 
Following this procedure inductively, we set f2 = Qoo which will look like figure 3 

The red area is transported to xi, disconnected and has infinitely many compo- 
nents. 
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3. The sphere 

In the case of M — S", the n-dimensional sphere, the picture changes rather 
drastically. It seems that there is just not enough space on the sphere to produce 
disconnected sets. We have 

Theorem 2. Let M = S"' and consider the optimal transport problem from ji to v 
with cost function Cp{x, y) = ^dP{x, y). Then, for any i with Aj > and any p>l 
the setT~^{xi) is simply connected. 

For the proof we need two lemmas. Let us fix some notation first. For z G 
we denote by z' its antipode. For z,y G S"' we denote by Cz,y the closed geodesic 
connecting z and y. Then we have 

Lemma 2. Let C be any circle in S", i.e. C = dB{q,r) for some r > and 
q € S". Then, for any i there exist two points on C, say xf and xj , such that 
the function (z) attains its maximum and minimum on C in xf respectively xj . 
Moreover, {x^,x{} = CC\Cxi,q. In particular, CdHj is either empty or connected 
for any choice of i and j . 

Proof. If $j is constant on C the statement is clear, as we can choose any two points 
on C. If $i is not constant on C, by continuity of the distance function, it is clear 
that there are two such points as claimed. The great circle Cx^.q divides C into two 
halves on which $j is symmetric. Thus, it suffices to show that <I>i is monotonely 
decreasing from x" to xj, otherwise said, that the function Di{z) = d{xi,z) is 
monotonely increasing from a;" to xj . However, this is a direct consequence of 
the spherical law of cosine for the geodesic triangle A(xi,q,z) for some z E C. 
Parametrize the points z on C by the angle relative to the geodesic from q to 
Xi. Then we have 

cos{d{z, Xi)) = cos(r) cos(d(a;i, q)) + sin(r) sin(d(a;i, q) cos{az). 

If we now walk on C from a;f to x{ the angle is increasing from to w. 
Therefore, Di{z) has to be increasing as well. The last statement of the lemma 
is now a direct consequence of the monotonicity of the two functions and $j 
between its respective extremal points. Indeed, parametrize C by S^. Then, on C 
both $i and $j are 27r-periodic and symmetric with respect to Xi respectively Xj 
and have the same shape. They are shifted versions of each others. Thus, proving 
the claim. □ 

Lemma 3. If x', G m then H] = 5*". 

Proof, x'l G i/j means that ^i{x'j) > ^j{x'i). We have to show that for arbitrary z G 
5" we have ^i{z) > ^j{z). For given z we consider the circle C = dB{xj , d{xj , zj) 
around Xj, where $j has constant value ^j{z). From the last lemma we know 
that there is a point y & C such that ^i{y) < ^i{w) for all w e C. Furthermore, 
y G C (1 Cxj,xi- Hence, it suffices to show that > ^j{y)- For this, it is 

sufficient to show that > ^j{q) for all q G Cx^.xi —■ C. Parametrizing C as 

we see that both $j and are 27r-periodic. Moreover, is just a shift of 
that is 

= ^i{q + d{xi, Xj)) + bj - bi, 
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where the sum inside the brackets is a sum in considered as a group and the 
sum outside the brackets is a usual sum in M. Thus, > ^j{q) for all q ^ C iff 

<^i{q) - + d{x^,Xj)) > bj - bi. 

By the special form of <i>i, the left hand side is bounded below by 

<^i{x'^) - ^^{x'^ + d{xi,Xj)) = <l>i{x'^) - <i>j{x',) - bj + bi> b, - bj > > bj - h, 

as 

b,-^ = ^,{x\)><^j{x'^>b,^^ 

implying bi > bj . 

□ 

This enables us to prove the theorem. 

Proof of theorem\^ Assume y,z ^ T^^{xi), i.e. ^i{y) > ^jiy) for all j ^ i and 
for z likewise. Consider the circle C going through y, z and x\. As iJj n C =: G is 
not empty (y, z G -ffj n C) it is connected and y and z have to lie in this connected 
segment of C. Thus, either x\ ^ G oy x\ ^ G. In the former case however we 
have ^i{x[) > ^j{x[). By the previous lemma, this means 7?] = S'", in particular 
C C Hj. Therefore, in either case the segment of C connecting z and y without 
passing through x'^ is transported to Xi . As this is independent of j this shows that 
T~^{xi) is pathwise connected and therefore connected. 

Assume that T~^{xi) is not simply connected. Then there with k ^ 

j ^ i ^ k such that any path joining T~^(xj) and T^^(xk) has to pass through 
T^^{xi). Without loss of generality we can assume that x'^ £ T^^{xj) (if there is no 
such j the whole sphere is transported to Xi). Take any z € T~^[xk) and consider 
Cz.x'.- By assumption, on either segment of Cz^x'. connecting z and x[ there are 
subsegments mapped to x^, say Gi and G2. From the lemma above we know that 
HlnCz^x' is connected. As Gi and G2 are mapped to Xi they have to be connected 
by a segment of C^^x'., say G3. Then either z £ G3 or a;^ £ G. The latter again 
implies that C^^x'. = G3 and, in particular, we thus have z G G3. But then z cannot 
be mapped to Xfc, a contradiction. □ 

4. General manifolds and "p — 1" 

We now want to consider the problem on a general complete compact Riemannian 
manifold. For the case p — I the situation is very nice for all manifolds and we get 
a good geometric description of the cells of the Laguerre tessellation. Due to the 
triangle inequality, we have the following 

Theorem 3. Let M be a compact connected complete Riemannian manifold without 
boundary. Consider the optimal transport problem from ^ to v with cost function 
ci{x, y) = d{x, y), the geodesic distance. Then, for any i with > the set T~^(xi) 
is simply connected and starlike with respect to Xi. 

Proof. We claim that the halfspace Hj is starlike with respect to Xi. Then, if 
z G T~^{xi) = C\j-ii Hj we have z € Hj for all j. As Hj is starlike with respect to 
Xi, the geodesic jt connecting z and Xi lies entirely inside H^. And this holds for 
all j. Thus, {74} C r\j^i^j ^^'^ T~^{xi) is not only connected but also starhke 
with respect to Xi. Hence, it suffices to prove the claim. 
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Assume T{xi) ^ xi. Then, there is a j ^ i such that T{xi) — xj, i.e. ^j{xi) = 
—d{xi,Xj) + bj > bi ~ ^i{xi). Then, we have for any q ^ M 

-d{q,Xj) + bj > -d{q, Xi) - d{xt,Xj) + bj > -d{q,Xi) + bi. 

This impHes, that T^^{xi) = contradicting the assumption of Xi > 0. Thus, 
T{xi) = Xi. 

Take any w G T^^(xi) (hence, ^i{w) > ^j{w) for all j ^ i) and q ^ M such 
that d(xi, w) = d{xi, q) + d{q, w), i.e. q lies on the minimizing geodesic from Xi to 
w. Then, we have for any j ^ i hy using the triangle inequality once more 

—d{q,Xi) + bi = —d{xi,w)+d{q,w) + bi 

> —d{xi, w) + bi + d(w, Xj) — d{q, Xj) 

> -d{q, Xj) + bj, 

which means that ^i{q) > ^j{q) for all j ^ i. Hence, q G Hj^i^j proving the 
claim. □ 

Remark 1. Note that in this case the transport map is unique even though we are 
working with the cost function ci{x,y) — d{x,y). The reason is that we transport 
the mass to a discrete measure. The sets dH^j = {x ^ M : = <I>j(a;)} are null 

sets. Thus, their (countable) union is a null set and for any z G M\\^^ ■ dHj there 
is a unique transport plan. 

However, as mentioned above, already for the case p = 2 on a compact subset 
of the hyperbolic space the cells T^^{xi) can be disconnected (see Example 3.8 in 
[KM07| ). Results of various simulations hint at this not being special for p — 2. 
We believe, that for the hyperbolic space, this holds for all p > 1. In the case of a 
compact manifold with positive and negative curvature, we can be even worse off. 

Example 2. Let p > I be given. Consider two points on the .sphere which are the 
antipodes of each other, say x and x' . Fix the transport map T = exp(V$) with 
<i> = max{$i, <i>2} and ^i{y) = -^dP{x, y) + bi and ^2{y) = -^dP{x',y), as above. 
Assume, that bi is such that $i(a;') = $2(2^'). In particular, we have T{S"') = x. 
We now want to deform the sphere into a manifold M with positive and negative 
curvature while fixing the transport map T ( the metric of the manifold will change 
and therefore we will end up with another distance function, but we will suppress 
that in the notation). 

Fix a geodesic connecting x and x' ( that is one half of a great circle connecting 
X and x' ). The x' -equator divides the geodesic into two segments. One of the 
segments, say G, has the property that all of its points are closer to x' than to x. 
Let us fix a point g G G not equal to x' . By adding a very thin spike to 5" around 
g and making it, if necessary, very long, we can achieve, that the top of the spike is 
transported to x' . The reason is, that the cost function c{z,y) — ^dP{z,y) is strictly 
convex and the points on the spike are closer to x' than to x (take k arbitrary, p > 1, 
then, there is R such that for all r > R we have —(r + l)P-\-K< —{r— lY ). Then, 
in a similar manner to example Q] we can add a sequence of spikes not touching 
each other, getting smaller and smaller and converging to a "zero-spike" at x' (see 
figure 4)- In- order that this yields infinitely many components of the cell T~^(x') 
we need the assumption <^i{x') = $2(3:'). 
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Figure 4. Example [2] 

Thus, a natural question is, what is the geometry of the cells T^^{xi) in an 
arbitrary compact connected Riemannian manifold with curvature bounded below 
by a strictly positive constant. Rough estimates on deformation of the sphere 
suggest that the cells should at least be connected. However, so far we were just 
able to prove connectedness in the case oi p = 2, curvature bounded below by zero 
under the additional assumption that T{xi) = Xi. 
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